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9C continued: results from a deeper radio-source survey at 15 GHz 
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ABSTRACT 

The 9C survey of radio sources with the Ryle telescope at 15.2 GHz was set up to survey 
the fields of the cosmic microwave background telescope, the Very Small Array. In our first 
paper we described three regions of the survey, constituting a total area of 520 deg^ to a 
completeness limit of « 25 mJy. Here we report on a series of deeper regions, amounting 
to an area of 115 deg^ complete to « 10 mJy and of 29 deg^ complete to w 5.5 mJy. We 
have investigated the source counts and the distributions of the 1.4 to 15.2 GHz spectral index 
(aP4^) for these deeper samples. The whole catalogue of 643 sources is available online. 

Down to our lower limit of 5.5 mJy we detect no evidence for any change in the differen- 
tial source count from the earlier fitted count above 25 mJy, «(5') = 51 (S'/Jy)^^ '^ Jy^'sr^^ 

We have matched both our new and earlier catalogues with the NRAO VLA Sky Survey 
(NVSS) catalogue at 1.4 GHz. For samples of sources selected at 15.2 GHz, in three flux 
density ranges, we detect a significant shift in the median value of a\^^^; samples with higher 
flux densities have higher proportions of sources with flat and rising spectra. We suggest 
that this observed shift is consistent with a model containing two distinct source populations 
having differently sloped source counts. Samples selected at 1.4 GHz contain significantly 
smaller proportions of sources with flat and rising spectra. Also, in our area complete to 
« 10 mJy, we find 5 sources between 10 to 15 mJy, amounting to 4.3 per cent of sources in 
this range, with no counterpart in the NVSS catalogue. These results illustrate the problems 
inherent in using a low frequency catalogue to characterise the source population at a much 
higher frequency and emphasise the value of our blind 15.2-GHz survey. 
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1 INTRODUCTION 

As explained in our first paper (Waldram et al. 2003), hereafter Pa- 
per I, the 9C survey with the Ryle Telescope (RT) at 15.2 GHz 
was designed specifically as part of the observing strategy of the 
Very Small Array (VSA) (Watson et al. 2003). Foreground radio 
sources are a major contaminant in the cosmic microwave back- 
ground (CMB) observations of the VSA at 34 GHz and the prime 
purpose of the 9C survey has been to provide a catalogue of radio 
sources in each of the VSA fields (Taylor et al. 2003; Cleary et al. 
2005, 2008). This has meant that the positions of the 9C fields and 
their depth in different areas have been determined by the VSA ob- 
serving programme. For the VSA compact array {(. = 150 to 800) 
we covered a total area of 520 deg^ to a completeness limit of 
25 mJy but for the extended array ((. = 300 to 1500) and now the 
super-extended array (i = 600 to 2500) deeper surveying over some 
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smaller regions has been required. It is these regions which are de- 
scribed here. 

In addition to its specific application to the VSA observations, 
the 9C survey has proved to be of much wider interest. There have 
been other blind surveys at high radio frequencies, notably with 
the Australia Telescope Compact Array (ATCA) at 1 8 and 20 GHz 
(Ricci et al. 2004; Sadler et al 2006; Massardi et al. 2008), but 9C 
was the first such survey of any extent and is the only one to reach 
such low flux densities. It provides an important opportunity for 
exploring the radio-source population at high frequencies. Samples 
of 9C sources have been followed up with simultaneous multifre- 
quency observations at 1.4, 4.8, 15.2, 22 and 43 GHz (Bolton et al. 
2004) and there have been two further papers based on these data: 
Bolton et al. 2006a, which describes results from 5-GHz MERLIN 
and VLBA observations, and Bolton et al. 2006b, which examines 
variability at 15.2 GHz. 

As a blind survey at 15.2 GHz, 9C is particularly valuable in 
the study of source counts at high radio frequencies. Point sources 
are a significant foreground in the CMB observations of many 
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Table 1. The areas complete to « 10 mjy 



Field 


RA 


(- 


^) to RA 


DecC 


") to Dec C 


"') 


0020+2947 


00 


06 


48.2 


00 


33 04.6 


+26 


37 


55 


+32 


55 


20 


0303+2629 


02 


49 


36.5 


03 


16 03.3 


+24 


26 


19 


+28 


32 


04 


0731+5427 


07 


22 


05.0 


07 


40 01.2 


+53 


09 


29 


+55 


44 


20 


0938+3218 


09 


29 


38.0 


09 


46 08.0 


+30 


43 


45 


+33 


51 


55 


1228+5301 


12 


17 


08.5 


12 


38 43.9 


+51 


27 


16 


+54 


33 


45 


1535+4305 


15 


22 


16.7 


15 


47 10.1 


+40 


43 


16 


+45 


25 


52 


1735+4215 


17 


29 


40.3 


17 


43 57.2 


+41 


13 


24 


+43 


17 


30 



instruments (e.g. the Planck Surveyor) and, even if the brighter 
sources can be identified and subtracted from the data, there re- 
mains the problem of residual confusion from the fainter, unsub- 
tracted source-population. This effect can be evaluated only from a 
knowledge of the source counts below the subtraction limit. Blind 
surveys at frequencies of 30 GHz and above are extremely prob- 
lematical because of the small field of view of the available tele- 
scopes, so one approach has been to identify sources from a lower 
frequency survey, with a measured source count, and follow them 
up at the higher frequency. 9C, at 15.2 GHz, is comparatively close 
to the target frequencies and, therefore, has a considerable advan- 
tage over lower-frequency surveys. 

We have used this method to make empirical estimates of the 
source counts at higher frequencies (including those in the lower 
Planck bands) by taking the spectral index distributions from the 
Bolton data together with our known 15.2-GHz source count (Wal- 
dram et al. 2007). Comparison with the theoretical model of de 
Zotti et al. (2005) shows good agreement below 43 GHz but at 
higher frequencies there is some divergence towards higher flux 
densities, indicating that our predictions imply fewer flat spectrum 
sources. However, the results of Sadler et al. (2008) who have made 
simultaneous observations with ATCA of a sample of sources at 20 
and 95 GHz, do appear to confirm our own estimates of the count 
at 90 GHz in the common flux-density range. 

More recently. Mason et al. (2009) have made a 31 -GHz sur- 
vey of the fields of the Cosmic Backgroimd Imager (CBI, Mason 
et al. 2(X)3) by following up known extragalactic sources from the 
NVSS survey at 1.4 GHz (Condon et al. 1998). They have esti- 
mated the spectral index distribution from 1 .4 to 3 1 GHz and hence 
the 31-GHz source count in the range 0.5 to 10 mJy, which has en- 
abled them to make an assessment of the contribution of unresolved 
point sources to the CBI power spectrum. Our results here provide 
a useful comparison with their work. 

Also, a subset of the sources in our present catalogue has now 
been followed up at 30 GHz with observations made using the 32- 
m, Torun telescope (Gawronski et al., in preparation). 

In this paper, we first explain the choice of fields and our meth- 
ods of observation and data analysis, together with estimates of 
completeness (Sections 2, 3 and 4). In Section 5 we present the 
extended source counts and in Section 6 investigate the correlation 
of our source lists with the NVSS survey. Section 7 describes the 
source catalogue. Section 8 discusses our conclusions and Section 
9 our plans for future work. 



2 THE SURVEY AREAS 

The seven areas of the survey presented here are centred at: 
0020+2947, 0303+2629, 0731+5427, 0938+3218, 1228+5301, 
1535+4305, 1735+4215 (RA (•> '») Dec. {° '), J2000.0), as indicated 
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Figure 1. The deeper 9C .survey areas: an equatorial plane projection with 
the N. pole at the centre. The Dechnation circles are at intervals of 10° and 
the Galactic plane is shown. 



Table 2. The areas complete to « 5.5 mJy 



Field 


RA('' -)toRA('i ^1 


Dec ( ' ") to Dec (" ' ") 


0020+2947 


00 1 1 55.7 


00 20 08.0 


+27 09 39.0 +29 16 04 




00 20 08.0 


00 22 45.5 


+28 43 43.0 +29 16 04 




00 25 00.8 


00 30 54.6 


+31 53 24.0 +32 24 07 


0303+2629 


02 55 27.2 


02 58 26.5 


+24 56 03.0 +25 27 56 




02 58 26.5 


03 01 23.7 


+25 27 56.0 +25 58 27 




02 58 21.7 03 04 22.8 


+26 29 55.0 +26 59 31 


0731+5427 


07 22 05.0 


07 31 04.6 


+53 09 29.0 +53 40 27 




07 26 33.2 


07 40 01.2 


+53 40 27.0 +55 13 28 




07 22 05.0 07 31 17.7 


+55 13 28.0 +55 44 20 


0938+3218 


09 29 38.0 


09 40 46.9 


+30 43 45.0 +32 49 53 




09 43 01.8 


09 46 08.0 


+31 14 44.0 +3146 05 




09 40 46.9 


09 43 35.0 


+32 17 15.0 +32 49 55 




09 32 08.6 


09 35 23.3 


+32 49 53.0 +33 21 31 




09 37 38.7 


09 40 46.9 


+32 49 53.0 +33 21 16 


1228+5301 


12 21 19.2 


12 25 44.9 


+51 58 14.0 +52 30 36 




12 25 44.9 


12 30 01.0 


+51 58 14.0 +54 02 56 




12 34 22.2 


12 38 48.0 


+51 58 21.0 +52 30 43 




12 20 54.5 


12 25 44.9 


+53 3151.0 +54 02 56 


1535+4305 


15 34 28.4 


15 44 07.5 


+40 43 16.0 +42 49 11 




15 44 07.5 


15 47 10.1 


+42 18 13.0 +42 49 11 




15 22 16.7 


15 28 44.7 


+42 48 44.0 +43 19 48 




15 22 16.7 


15 25 32.3 


+43 19 48.0 +43 51 07 




15 40 59.3 


15 44 46.9 


+44 23 59.0 +44 55 40 




15 24 45.0 


15 31 50.2 


+44 54 51.0 +45 25 52 


1735+4215 


17 29 40.3 


17 43 57.2 


+41 13 24.0 +42 15 51 




17 29 40.3 


17 33 16.0 


+42 15 51.0 +42 46 40 




17 36 50.1 


17 43 57.2 


+42 15 51.0 +42 46 40 
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Figure 3. A typical sensitivity map corresponding to a 9 x 8 hexagonal ar- 
ray of pointing directions. White indicates high sensitivity. The coordinates 
are in map pixels, where one pixel represents 8x8 arcsec- . 

in Figure[T] They are situated away from the Galactic plane and are 
widely spaced in RA, their positions having been determined by the 
choice of fields for the observations of the VSA. The RA and Dec. 
ranges are given in Table[T] The total area amounts to 114.7 deg^ 
and within this our source list is complete to ~ 10 mjy. In each 
field there are also some much more sensitive areas; this was partly 
because deeper surveying was required near the centre of the VSA 
primary beam and partly because on some days there were par- 
ticularly favourable observing conditions. For the purpose of this 
paper, in order assemble a complete sample of a useful size, we 
have selected a number of sub-areas with ~ 5.5 mJy completeness, 
which form a total area of 29.1 deg^ (Figure[2l(. These are of vari- 
ous shapes and sizes, so in Table[2]we have described them in terms 
of small constituent areas bounded by specific RA and Dec. ranges. 
Although our completeness limit in these regions is ~ 5.5 mJy, we 
have detected many fainter sources, the faintest being only 1 mJy. 
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Figure 4. Plot of the ratio (Sr/5p) of raster flux density to pointed flux 
density versus Sp. for the deeper areas. The estimated completeness limit is 
5.5mjy. (Where the ratio is shown as zero, the corresponding source lies 
in an area of a raster map which is hard to interpret. A pointed observation 
has been made but there is no definitive raster flux density.) 



The source extraction was also implemented much as before. 
Possible sources were identified as peaks from interpolation of the 
raster map and then each of these was followed up with a short 
pointed observation of 10-15 minutes, either to establish a reliable 
flux density or to eliminate it as a false detection. The only change 
in our method was that, while we continued to take 3cj as the pixel 
cut-off on the raster map, we now took 4.5CT, rather than 5.0<J, as 
the cut-off for the interpolated peak, in order to find as many of the 
weaker sources as possible. 



3 OBSERVATIONS AND DATA ANALYSIS 

For our observations we have continued to use a rastering tech- 
nique similar to that described in Paper I. This is necessitated by 
the relatively small area within the FWHM of the RT primary beam 
(0.01 deg^) compared with the areas of the VSA fields to be sur- 
veyed. As before, we used the RT's E-W array of five aerials to 
give a resolution of 25 x 25cosec8 arcsec^ and set up a hexago- 
nal pattern of pointing directions with a spacing of 5 arcmin. But, 
in order to reach the required sensitivity, we now scanned only 72 
pointings during a 12h observation, covering an area of sky of only 
~ 0.3 deg^. We also needed to integrate several observations of the 
same area, sometimes as many as six for the deeper regions, if the 
weather were unfavourable. Separate so-called 'constituent' maps 
were made from the accumulated data from the individual point- 
ings and, in an improvement on our earlier method, these were then 
CLEANed before assembling them to form the 'raster' map. The 
method for combining the constituent maps was the same as before 
and produced a similar variation in sensitivity (of ~ 20 %) over the 
corresponding raster map, as illustrated in Figure [3] We made an 
estimate of the mean noise, a, over each raster map and continued 
to integrate data until we had reached the required sensitivity. Over 
all the areas we reached o < 2 mJy and in the deeper areas o < 1 
mJy, though there was considerable variation within these. 



4 COMPLETENESS 

Since the original candidate sources were found from raster maps 
with a range of noise levels, assessment of the completeness limits 
of our source lists is not straightforward. Although we can be sure 
that every source in our catalogue is genuine, since it has been fol- 
lowed up with a pointed observation, it is difficult to estimate the 
number of undetected sources at any flux density level. We have 
approached the problem in the following way. 

Taking the deeper areas first, we have plotted the ratio, 
(5r/5p), of the raster flux density to the pointed flux density ver- 
sus 5p, as in Figure |4] We find that, for the brighter sources, in 
the range > 10 mJy where we can assume completeness, there is a 
scatter about a median ratio of ~ 0.89. This off-set of ~ 10 per cent 
in the median value is probably attributable to pointing problems; 
in particular, there was a necessary compromise between keeping 
as many of the observed data samples as possible and discarding 
those with excessive pointing errors. For the weaker sources the ra- 
tio rises progressively with decreasing flux density because at these 
levels the peaks detected on the raster maps are preferentially those 
boosted by a positive local noise contribution. In order to estimate 
the completeness limit, we have investigated the change in the rel- 
ative numbers of sources above and below the median line as the 
pointed flux density decreases from 10 mJy. We find that there is 
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Table 3. The data used in the combined count 



bin start 


bin end 


number 


area 


S/Jy 


5/Jy 


N 


sr 


0.0055 


0.0063 


23 


0.008856 


0.0063 


0.0078 


21 


0.008856 


0.0078 


0.0100 


21 


0.008856 


0.0100 


0.0110 


38 


0.034939 


0.0110 


0.0129 


43 


0.034939 


0.0129 


0.0152 


40 


0.034939 


0.0152 


0.0198 


40 


0.034939 


0.0198 


0.0250 


39 


0.034939 


0.025 


0.030 


80 


0.15841 


0.030 


0.040 


104 


0.15841 


0.040 


0.060 


103 


0.15841 


0.060 


0.100 


94 


0.15841 


0.100 


0.200 


49 


0.15841 


0.200 


0.500 


27 


0.15841 


0.500 


1.000 


8 


0.15841 




S(Jy) 

Figure 6. The source count for the sources in the deeper areas, showing the 
lower completeness limit of r; 5.5mJy and the shortfall above 100 mjy 



fitted count from Paper I 



S(Jy) 

Figure 5. The source count for the sources in the total areas, showing the 
lower completeness limit of « lOmJy and the shortfall above 100 mJy 



520 deg' 
115 
29 

fitted count from Paper I 



S(Jy) 

Figure 7. Differential count extended to 5.5 mJy 



no significant difference in the numbers down to 5.5 mJy but be- 
low this level there are significantly more sources with ratios above 
0.89 than below. We therefore assume we have completeness to ap- 
proximately 5.5 mJy, which means there are 135 sources above the 
limit and 8 1 below. 

For the total areas, we apply similar arguments using higher 
flux densities and conclude we have a completeness limit of ap- 
proximately 10 mJy, giving us 307 sources above the limit and 336 
below. 

We should also explain that within the areas presented in this 
paper there is a shortfall in the number of sources above 100 mJy 
as compared with the number expected from the source count in 
Paper I (see next section). This is because these VSA fields were 
deliberately selected to contain as few bright sources as possible. 



5 SOURCE COUNTS 

The source count for the sources in the total areas is shown in Fig- 
ure |5] which indicates the lower completeness limit of 10 mJy and 
also the shortfall above 100 mJy mentioned in the previous section. 
Figure|6]is a similar plot for the sub-set of these sources which lie 
in the deeper areas. 

Using these counts, we are now able to extend the original 9C 
count to a lower flux density level. We have combined data from 
Paper I in the range 5 > 25 mJy with data from this paper in the 
range 10 < 5 < 25 mJy from the total areas and 5.5 < 5 < 10 mJy 
from the deeper areas (see Table |3j. As explained in the earlier 
paper, the effect of the random errors in the flux densities is negli- 
gible compared with the Poisson errors on N, but we have applied 
a correction for the bin widths in the same way as before. At the 
completeness limit of 5.5 mJy the survey continues to be limited 
by noise rather than confusion; we estimate that above this level 
there are typically 2000 synthesized beam areas per source. 
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Figure 8. Normalised count compai'ed with tlie de Zotti model 

Figure [7] shows the differential count ;2(5) and we see that, 
down to 5.5 mJy, there is no evidence for any significant change 
from the fitted count n{S) calculated in Paper I: i.e. 

In Figure [S] we compare our normalized count S^'^n{S) with that 
derived from the model of de Zotti et al. (2005) and note that our 
count appears to lie somewhat below their prediction in the lower 
flux density range. 

These two figures effectively update the plots in Figures 6 and 
7 of Paper I. We should mention here that we found retrospectively 
that the positions of two points in each of those figures are in error. 
These are the two open circles corresponding to the higher flux 
bins of the deeper survey. In fact, since they lie in the flux density 
range of the main survey, they do not affect any of the arguments 
or conclusions in that paper and can simply be ignored. 



6 CORRELATION WITH NVSS AT 1.4 GHZ 

In this section we present the results of our investigations but defer 
any discussion of their implications to Section 8. 



6.1 Matching the catalogues 

We have attempted to match the 643 sources in our catalogue with 
the catalogue from the NVSS 1.4-GHz survey, which has a resolu- 
tion of 45 arcsec (not dissimilar to our own) and a completeness of 
50 percent at ~ 2.5 mJy rising rapidly to 99 percent at 3.4 mJy. In 
the first instance, we searched for a counterpart within 40 arcsec of 
each of the 9C sources and then in every case inspected the corre- 
sponding NVSS contour plot, paying particular attention to sources 
which had multiple matches within 40 arcsec and to those where 
there was no match within 20 arcsec. For about 90 percent of the 
sources the matching was straightforward; the remainder, however, 
required further attention in order to make the appropriate correla- 
tion of flux densities. For example, some binary sources had a sin- 
gle entry in the 9C catalogue, but appeared separately in the NVSS 



0.1 r 



0.001 




0.01 0.1 1 

1 .4 GHz flux density (NVSS or FIRST) (Jy) 



Figure 9. Plot of RT 15.2 GHz pointed flux densities versus NVSS (or 
FIRST) 1.4 GHz flux densities, with a line indicating zero spectral index. 
The horizontal line is the estimated completeness limit. 



Table 4. Some statistics for the the spectral index distributions otj 



15.2-GHz flux 

density 
range (mJy) 


Number 

of 
sources 


Number of 
sources 
with 

01^4- < 


Number of 
sources 
with 
< 0.5 


"1.4 


5.5 <5<25 


111 


11 (10 %) 


40 (36 %) 


0.79 


25 < 5 < 100 


381 


62 (16 %) 


190 (50 %) 


0.51 


S> 100 


84 


26 (31 %) 


60 (71 %) 


0.23 



catalogue and, conversely, there were some single NVSS sources 
which appeared as separate components in 9C. Where possible, we 
made use of the deeper FIRST 1.4-GHz survey (Becker, White & 
Helfand 1995), with its resolution of 5 arcsec, to assist in resolving 
ambiguities. 

Our correlation is illustrated in Figure |9] In all, out of the 643 
sources in our list, there were 21 with no counterpart in the NVSS 
catalogue, though two of these did appear in FIRST. Also, a num- 
ber of them were found to be coincident with one or two contours 
on the NVSS plots but clearly these were not significant enough de- 
tections to justify their inclusion as sources in the NVSS catalogue. 
(See sub-section 6.3 for further details.) 



6.2 Spectral index distributions: samples selected at 
15.2 GHz 

A histogram of the spectral index distribution aj^^^ for the 111 
sources with 5.5 < 5 < 25 mJy and lying in our deep areas is shown 
in Figure[TOl We define a by S oc v^". The figure also shows the 
spectral index distributions, over the same range of frequencies, for 
sources with 25 < 5 < 100 mJy and S > 100 miy; these distribu- 
tions contain 381 and 84 sources, respectively, and were originally 
included in Paper I. They are derived from an area of 520 deg^ com- 
plete to ~ 25 mJy. In Table|4]we have summarised some important 
features of the distributions; it is clear from this table and FigurefTOl 
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Table 5. Results from a Monte Carlo simulation, of one million realisations, to assess the significance of the shift in 
median spectral index with flux density. 



15-GHz flux 

density 
range (mJy) 


Number 
of 

simulated 
sources 


Min. 

"l.4 


Max. 
"l.4 


Mean 
a|4" 


Number of 
realisations 
with a [4^ 


Probability of 
drawing median 
more extreme 
than observed 


5.5 <5< 25 


ill 


0.166 


0.817 


0.477 


> 0.79 is 14 


2.8 X 10-5 


25 < 5 < 100 


381 


0.311 


0.685 


0.474 


>0.51 is 212833 


0.426 


S> 100 


84 


0.081 


0.879 


0.478 


< 0.23 is 161 


3.2 X ur"^ 



-0.5 □ 

Spectral Index 



-0.5 0.5 
Spectral Index 



S 15 ^ 



-0.5 0.5 

Spectral Index 



-0.5 
Spectral Index 



Figure 10. The otl'^' spectral index distributions for samples selected at 
15.2 GHz in three flux density ranges. From the top: 5.5 < 5 < 25 mJy, 
25 < S < 100 mJy and S > 100 mJy, with median values of 0.79, 0.51 and 
0.23 respectively. We define a by 5 « v"". (The dotted box in the top plot 
indicates a limiting value of a for a source not in NVSS) 



Figure 11. The aj^^- spectral index distributions in Figure [Tol combined 
into one distribution. 5 « V"" 



that, as one moves to lower flux densities, the proportion of sources 
with flat and rising spectra decreases. 

In order to assess the significance of the shift in median spec- 
tral index with flux density, we have combined the spectral indices 
of the 576 sources to form a single spectral index distribution, 
which we have plotted as a histogram in Figure [TT] We assume 
that there is an underlying spectral index distribution which has no 
dependence on flux density and that the combined data represent 
our best approximation to this distribution. For each of our flux 
density bins, we calculate the probability of drawing the observed 
number of sources from the combined distribution and finding a 
median spectral index as extreme as those observed, using a Monte 
Carlo simulation of 10^ realisations. For example, for the lowest 
flux density bin, for each realisation, we draw, at random, ill val- 
ues from the combined spectral index distribution and find the me- 
dian spectral index. We find 14 realisations, out of 10^, in which 
the this value is greater than or equal to the observed value of 0.79 
and conclude that the probability of drawing ill sources from the 
combined distribution and observing a median as or more extreme 
than 0.79 is 2.8 x 10"^ 

Table |5] provides a summary of the results from the simu- 
lations. It shows that the median spectral index observed for the 
middle flux density bin is not significantly different from that of 
the combined distribution of Figure [TT] However, the probabili- 
ties of having obtained median spectral indices more extreme than 
observed, for the two outer flux density bins, by drawing sources 
randomly from the combined distribution are very small. Conse- 
quently, we can reject our null hypothesis, that the underlying spec- 
tral index distribution has no dependence on flux density, with a 
high level of confidence. In other words, the shift in the median 
spectral index with flux density is significant. 
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Table 7. 9C sources without corresponding matclies in NVSS. Flux densi- 
ties for those sources appearing in the FIRST survey are provided, along 
with the spectral index based upon the FIRST flux density. For sources not 
appearing in FIRST, an upper limit to the spectral index is provided by as- 
suming a flux density of 3.4 mJy (the estimated 99 per cent completeness 
limit of NVSS) at 1.4 GHz. 



Figure 12. Histogram of the spectral index distribution a[ for a sam- 
ple selected from NVSS with 5i 4 > 32.95 mJy, showing the percentage of 
sources with flat and rising spectra. 5 « v^" 



6.3 Spectral index distributions: samples selected from 
NVSS at 1.4 GHz 

It is not possible to calculate a full spectral index distribution for 
a sample of NVSS sources complete to a level of 3.4 mJy, since 
many of the steep spectrum sources will not be detectable in our 
survey. We can, however, investigate the proportion of flat and ris- 
ing spectrum sources in the NVSS catalogue for different cut-offs 
in 1.4 GHz flux density, where we use 'flat and rising' to refer to 
sources with a^^^ < 0.5. 

We have selected all sources from the NVSS catalogue, ly- 
ing in the areas complete to 10 mJy at 15.2 GHz, with 5i,4 > 
32.95 mJy; in total, 644 sources are selected. These selection cri- 
teria ensure that any sources with ap^^ < 0.5 will have 5i5 2 > 
10 mJy, our completeness limit at 15.2 GHz. Of the 644 sources 
selected, 63 are found to have o}^^ < 0.5. In other words, we find 
that 9.8 per cent of sources with 5i 4 > 32.95 mJy have a^^^ < 0.5. 
We have plotted this flat and rising part of the spectral index distri- 
bution aj^4^ at 1.4 GHz as a histogram in Figure[T2l We have also 
taken a lower flux density cut at 1.4 GHz of 10 mJy (again in our 
total areas), which enables us to investigate sources with ap^^ < 
for a somewhat deeper sample. 

Similarly, we have taken appropriate NVSS flux density cuts 
in the deeper areas, complete to 5.5 mJy at 15.2 GHz, to enable us 
to investigate the rising, and flat and rising side of the distribution 
at lower flux densities. Our results are summarised in Table|6l 



6.4 Sources not in NVSS 

Table[7]liststhe21 sources in our total source list which do not ap- 
pear in the NVSS catalogue. Only one of these lies in the deeper 
areas (see the dotted box in the top plot of Figure llOl l and the 
others are all in the total areas. By dividing these 21 sources into 
three flux-density bins, as shown in Table [S] we can gain an in- 
dication of the proportion of 9C sources, with no counterparts in 
NVSS, as a function of flux density. We are complete only in the 
highest bin (10 < 5 < 15 mJy), for which we find 4.3 per cent of 
sources with no NVSS counterpart. Assuming that the proportion 
of sources with no counterpart increases with decreasing flux den- 
sity, our estimates of the proportions of 'missing' sources are likely 
to be under-estimates for the two lower flux density bins. This is 
because these bins are incomplete; sources with higher fluxes are 
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17:42:52.2 


-1-43:09:51 
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00:15:58.0 


-1-27:18:54 


3.1 




0.04 


00:31:12.7 


-1-27:12:58 
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0.05 


07:27:18.0 
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3.0 




0.05 
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-1-44:16:22 


2.6 




0.11 


Table 8. 9C sources with no NVSS counterpart as a function of flux density. 


Flux 


Number 


Number of 






density 


of9C 


sources with no 






range (mJy) 


sources 


NVSS counterpart 




S<4.2 


62 


8 (12.9 %) 






4.2 <5< 10 


281 


8 (2.8 %) 






10<5< 15 


117 


5 (4.3 %) 







preferentially detected in these incomplete bins, yet it is the brighter 
sources that are more likely to have an NVSS counterpart. 

We should note that several of the sources have been observed 
at other frequencies. For example, a number are close to sources 
in the 2MASS (Skrutskie et al. 2006) and ROSAT all-sky survey 
(Voges et al. 1999) catalogues. However, here we are concerned 
with the number of sources at 15 GHz that would be missed by 
reliance on the NVSS catalogue. 



7 THE SOURCE CATALOGUE 

A short section of the catalogue is shown in Table|9] The whole list 
is available online at http://www.mrao.cam.ac.uk/surveysltogether 
with the relevant tables: i.e. Tables [T] and |2] from this paper We 
have included all 643 sources, which means that a large number of 
them are below our completeness limits of ~ 10 mJy for the total 
areas and ~ 5.5 mJy for the deeper areas, the faintest being only 
1 mJy. Sources which lie within the deeper areas are marked D in 
the catalogue. In the comment column 'b' denotes a binary source 
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Table 6. Proportions of rising, and flat and rising sources in tlie NVSS catalogue. 



Flux density 

cut at 
1.4GHz(mJy) 


Number of 
NVSS 
sources 


9C Area 


Select 
sources 
with 


Number of 
sources 
selected 


32.95 
10.0 

18.12 
5.5 


644 
1891 
295 
770 


110 deg^ complete to 10 mjy 
110 deg^ complete to 10 mJy 
29 deg^ complete to 5.5 mJy 
29 deg^ complete to 5.5 mJy 


a|" < 0.5 
a p/ < 

aP4- < 0.5 
a 'P/ < 


63 (9.8 %) 
20(1.1 %) 
27 (9.2 %) 
11 (1.4%) 



Table 9. A section from the source catalogue. The letter 'D' means that the source lies in one of the deeper areas listed in Table|2] In the comment column 'b' 
denotes a binary source and 'e' an extended source. 



Source name 




RA J2000 


Dec. J2000 


Flux density 


Date 


Comment 








li m s 


o / // 


(Jy) 


yymmdd 




9CJ00 15+3 123 




00 


15:47.5 


31:23:35 


0.0127 


020505 


e 


9CJ00 15+27 18 


D 


00 


15:58.0 


27:18:54 


0.0031 


031124 




9CJ00 16+2804 


D 


00 


16:02.4 


28:04:28 


0.0027 


031124 




9CJ0016+3158 




00 


16:07.4 


31:58:53 


0.0092 


041016 




9CJ00 16+3239 




00 


16:11.0 


32:39:26 


0.0230 


041228 




9CJ00 16+2945 




00 


16:13.0 


29:45:11 


0.0156 


010815 




9CJ0016+3238 




00 


16:13.6 


32:38:37 


0.0264 


041228 




9CJ00 16+3 139 




00 


16:40.6 


31:39:04 


0.0110 


030203 


b 


9CJ00 17+3209 




00 


17:02.2 


32:09:19 


0.0345 


041016 




9CJ0017+2733 


D 


00 


17:04.6 


27:33:43 


0.0062 


031124 





and 'e' a source which appears extended relative to the synthesized 
beam size. 

We have not assigned individual eiror estimates to the source 
parameters. The position of a source is derived from a raster map, 
unless the follow-up pointed observation indicates a significantly 
different position, in which case that is substituted. However, since 
there are large overlaps between adjacent raster maps and a num- 
ber of sources appear on more than one map, it has been possi- 
ble to make a general estimate of the position accuracy. Taking the 
0020+2947 field, we find 39 sources with repeat position measure- 
ments. Apart from one extended source, in no case do the measure- 
ments differ by more than 10 arcsec and the median difference is 
~ 4 arcsec. The flux densities are in the range 4.6 mJy to 68.3 mJy. 
For the brighter sources we can refer to Paper I, where we matched 
a sample of 17 sources, all > 60 mJy at 15.2 GHz, with their coun- 
terparts in the Jodrell VLA calibrator survey (Wilkinson et al. 1998, 
and references therein), finding an accuracy for these sources of 
better than 3 arcsec. 

Our flux densities are derived from the pointed observations 
and the uncertainty in these is dominated by the random uncer- 
tainty in the flux calibration, which we estimate to be ~ 5 per cent. 
In the case of an extended source, an estimate is made of the inte- 
grated flux density. Since many of the sources are highly variable, 
the dates of the pointed observations are also included in the cata- 
logue. 



8 CONCLUSIONS 

By developing our techniques of observation and analysis, we have 
succeeded in extending the 9C survey to a completeness limit of 
~ 10 mJy in an area of 115 deg^ and of ~ 5.5mJy in an area 
29 deg^. As a blind survey at 15.2 GHz to this depth it provides 
a unique resource for the study of the extragalactic source popula- 
tion at high radio frequency. It is particularly valuable for the in- 
vestigation of the foreground point-source contamination in CMB 
observations, not only for the VSA but for a wide range of other 
CMB experiments. 



8.1 Source Counts 

In extending our differential source count to 5.5 mJy we detect no 
evidence of a change from the fitted power law calculated in Pa- 
per I: n{S) = 51(5/Jy)^^ '^ Jy^'sr^'. We have compared our nor- 
malized count S^'^n{S) with that derived from the model of de Zotti 
et al. (2005) and find that our count appears to lie significantly be- 
low their prediction in the lower flux density range. 



8.2 Spectral index distributions 

In discussing the a\^^^ spectral indices, we distinguish between the 
properties of a flux-limited sample selected at 15.2 GHz and one 
selected at 1.4 GHz. The spectral index distributions are necessar- 
ily different and the one cannot be inferred from the other in any 
simple way. 



In addition to the catalogues presented in this paper and in 
Paper I, there are source lists from further areas of the 9C survey 
yet to be published. We plan to make these available online in due 
course. 



8.2.1 1 5. 2-GHz samples 

We find that for samples selected at 15.2 GHz in three flux-density 
bins there is a significant decrease in the median a^^^ and increase 
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in the proportion of flat and rising sources with increasing flux den- 
sity. We cannot identify any observational bias that might produce 
this effect, but we can compare our results with those of Bolton et 
al. (2004) who made simultaneous multifrequency observations at 
1.4, 4.8, 15.2, 22 and 43 GHz of flux-limited samples of sources se- 
lected from 9C at 15.2 GHz. They found that raising the flux limit 
of the 9C sample produced a significant increase in the fraction of 
rising-spectrum sources. Their analysis differed from ours in that 
each of their samples had only a lower flux density limit and no 
upper limit and also they classified the sources in terms of their 
1.4 to 4.8 GHz rather than their 1.4 to 15.2 GHz spectral indices. 
However, they also showed that the 1.4 to 4.8 GHz spectral indices 
from simultaneous measurements were highly correlated with the 

1.4 to 15.2 GHz spectral indices obtained by simply matching 9C 
with NVSS, which is what we have done here. We believe we have 
observed a significant correlation between spectral index distribu- 
tion and flux density and that this does appear to be consistent with 
the work of Bolton et al. (2004). 

In order to assess whether such a change in spectral index dis- 
tribution is physically reasonable over such a relatively small range 
of flux density, we have constructed the following very simplified 
model. We postulate two different populations of radio sources, 
one flat spectrum and one steep spectrum. We model these two 
populations as having differential source counts ^ °= with 
P flat = —1.7 and ^sieep = —2.5. The number counts are normalised 
such that we have equal numbers of sources at 40 mjy. We gener- 
ate source counts between 5.5 and 500 mJy as shown in Figure [73] 
the source count from the combined distribution is consistent with 
that found from our data. We then assume that the the flat spectrum 
sources have spectral indices drawn from a Gaussian distribution 
with mean Ufiai = 0.2 and a{afia,) = 0.3; and that the steep spec- 
trum sources have spectral indices drawn from a Gaussian distribu- 
tion with mean dsteep = 0.8 and oio-steep) = 0.3. Using this model 
we generate the expected spectral index distributions for flux bins 

5.5 < 5 < 25 mJy, 25 < 5 < 100 mJy and 5 > 100 mJy, shown in 
Figure [T4l 

We stress that we have not attempted to optimise this highly 
simplified model in order to fit our data - a full physical model 
of the sort proposed by de Zotti et al. (2005) contains more than 
two populations and uses a more sophisticated parameterisation of 
the individual source counts than a simple power law. This model 
does, however, demonstrate that it is possible to produce a signif- 
icant shift in the spectral index distribution over a range of two 
orders of magnitude in flux density due to real astronomical causes 
rather than being generated by either systematic errors or selection 
effects. We expect that the spectral index distribution of a sample 
of sources selected at low frequency, e.g. 1.4 GHz, will show a 
much smaller variation with flux density owing to the smaller pro- 
portion of flat spectrum sources in such a sample. This expectation 
appears consistent with the spectral index distribution between 1.4 
and 30 GHz, for a sample selected at 1.4 GHz, observed by Mason 
et al. (2009). 

8.2.2 1.4-GHz samples 

As we have seen, it is not possible to calculate a full spectral in- 
dex distribution for a sample of sources selected from NVSS at 

1.4 GHz but it is useful to explore the flat and rising side of the 
distribution. We find that, for complete samples of sources above 
32.95 and 18.12 mJy at 1.4 GHz, there are 9.8 and 9.2 per cent re- 
spectively with a\^^ < 0.5, and for complete samples above 10 and 

5.5 mJy at 1.4 GHz, there are 1.1 and 1.4 per cent respectively with 
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Figure 13. Simulated source counts generated from a model of two pop- 
ulations of radio sources which follow power law source counts. The dot- 
dashed line corresponds to flat spectrum sources; the dashed line to steep 
spectrum sources; and the solid line to the total source population for this 
model. 
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Spectral index 

Figure 14. Simulated spectral index distribution generated from a model of 
two populations of radio sources. The solid line shows the spectral index 
distribution for sources in the flux density range 5.5 < S < 25 mJy; the 
dashed line for sources in the flux density range 25 < 5 < 100 mJy; the dot- 
dashed line for sources 5 > 100 mJy. These simulations correspond to the 
results shown in Figure flol 



a 1^4^ < 0. We note that these values are very close to those found 
by Mason et al. (2009) for their 1.4 to 31 GHz spectral indices. 
We see that the proportions of sources with ap^^ < and 
< 0.5 in samples selected at 15.2 GHz (Table |4]l are sig- 
nificantly higher than those in samples selected at 1.4 GHz (Ta- 
ble |6j. This effect, of selecting higher proportions of flat-and- 
rising-spectrum sources the higher the frequency, has been recog- 
nised before (see, for example. Peacock & Wall, 1981) and is as- 
cribed to a bias against detecting the increasing number of steep- 
spectrum sources whose spectra have peaked below the selection 
frequency and which therefore tend to fall below the flux-density 
cut-off. 
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8.3 Sources not in NVSS 

An important result from our matching with NVSS is the number 
of sources that we have detected at 15.2 GHz which do not appear 
in the NVSS catalogue, and, in particular, that in the flux density 
range of 10 to 15 mjy there are as many as 4.3 per cent not in 
NVSS. This is especially relevant to work that relies upon use of 
the NVSS catalogue to identify sources at high frequency. It is very 
likely to be significant both in the estimation of high frequency 
source counts and in the selection of individual sources for pro- 
jecting out of high frequency CMB data. Our work indicates that 
the proportion of sources that could be missed by reliance upon 
the NVSS catalogue might be as high as 4.3 per cent at low flux 
densities. 



9 FUTURE WORK 

We see that all our results are relevant to the estimation of extra- 
galactic source populations at high radio frequencies and of the 
consequent effects of foreground sources on CMB observations. 
In particular, there is the problem of the so-called 'CBI excess' 
(Mason et al. 2003), which is the apparent existence of an excess 
of power over intrinsic CMB anisotropy on small angular scales 
(/ > 2000). The question of whether this can be ascribed to con- 
tamination from the point-source foreground is discussed fully in 
Mason et al. (2009) and Sievers et al. (2009). 

We intend to collaborate with Mason et al. to explore further 
the comparison of the results in their paper (2009) with ours here 
and, if possible, to follow up our deeper source list with simul- 
taneous observations at 15 and 30 GHz. Since 9C is so close in 
frequency this should provide more useful information on the deep 
source counts at 30 GHz which are critical for interpreting the CBI 
excess mentioned above. 
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